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Abstract 



A transition in a spheroidal particle from the paraelectric to the ferroelectric 
phase as well as dynamic susceptibility are studied without approximation in 
the paraphase. It is assumed that the surface charge is compensated and the 
boundary condition for the polarisation is P — 0, i. e. with zero polarisation 
at the surface of the particle. There is an infinite number of resonance frequen- 
cies in the dynamic dielectric function within the quasistatic approximation. 
The paraphase properties of the dielectric response of the particle are discussed. 
The transition temperature decreases with decreasing diameter d of the particles 
inverse quadratically. There exists such a critical diameter that for the parti- 
cles with the diameter below the critical one the ferroelectric phase is absent. 
Comparison of the experiment with theoretical results is carried out. Introduc- 
tion of a dead layer thickness leads to a very good agreement of theory with 
the experiment for PbTiO^ and to a good agreement with the experimentfor 
BaTi0 3 . 



2 



1 Introduction 



Disordered ferroelectrics and antiferroelectrics have infrared, microwave and 
Raman spectra with their characteristic features which are not present in the 
ordered ferroelectrics and antiferroelectrics. Disorder in these materials may be 
of a different type. Stoichiometric order-disorder type systems are the best stud- 
ied materials of the mentioned type [I], [2] and pj]. Other type of order-disorder 
system are those materials in which two types of ions interchange their posi- 
tions at high temperatures. At lower temperatures there exists a freezing of their 
positions. When the freezing is theoretically very slow, this may be achieved 
practically by an annealing process at which the ions order, see in [2] and the 
references therein. The low-temperature phase is the ferroelectric phase. A sim- 
ilar type of materials are crystal-glass microcomposites. The nonstochiometric 
solid solutions in which one component is a ferroelectric type material and the 
other is a non-ferroelectric, or antiferroelectric type material form another class 
of order-disorder systems with ferroelectric or antiferroelectric properties. Their 
structural form or the form of clusters with different dielectric properties are 
usually of the nanocomposite form. New methods of the preparations of mate- 
rials have lead in the last years to a possibility to obtain nanostructural type of 
the order-disorder materials [2]. 

We studied the dielectric response of microcomposites of the ferroelectric- 
dielectric type in our papers [I] and [5]. A two-phase composite of ferroelectric- 
dielectric particles was studied as concerning the dielectric response within the 
quasistatic approximation. Maxwell-Garnett theory, effective medium theory 
and Bergman representation enabled us to calculate the dielectric function for 
a microcomposite. The dynamic and static properties of such composites have 
been reported in these papers. There is present a distribution of polar modes 
due to their presence in the original ferroelectric component in the bulk form and 
due to the presence of geometrical resonances. New low-frequency peaks appear 
due to the phenomenon of the geometric resonance. Their properties are de- 
pendent on the distribution of ferroelectric and non-ferroelectric clusters which 
percolate the material. The soft mode below the percolation transition within a 
state where at least one infinite ferroelectric cluster is present does not change 
its frequency depending on the concentration of the non-ferroelectric material 
whatever its strength changes. The soft mode disappears above the percola- 
tion transition, where there is no infinite ferroelectric cluster present. Thus 
the soft mode becomes a hard mode. Let us note that presence of the mul- 
tirelaxation phenomenon in the ferroelectric-dielectric composites of the type 
mentioned above has its origin in the presence of ferroelectric clusters with a 
distribution of their geometrical and topological properties. In materials of the 
order-disorder type with a regular crystal structure there is a possibility to find 
multirelaxation phenomenon due to other reasons [6], [7j and [8j where the com- 
plex dielectric function was calculated using semiempirical Bloch equations. A 
perturbation and a non-perturbation approaches were used to calculate a dielec- 
tric function, and the multirelaxation phenomenon was found. Its origin is due 
to presence of dispersion-less excitations in ferroelectric order-disorder materials 
with an incommensurate phase present at a given temperature. Clusters with an 
incommensurate modulation of the ferroelectric order parameter were observed 
[9J not only in dipolar glasses mentioned in [3], but also in Zr , 9S Hf . 02TW4 
ceramics exhibiting an incommensurate-commensurate phase transition due to 
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occurence of small polar-like regions (the domains within the incommensurate 
structure surrounded by discommensurations which appear upon cooling down 
from higher temperatures below some critical one). It is interesting that there 
exists a similarity in response of magnetic quantum systems with incommen- 
surably modulated phases and order-disorder systems with an incommensurate 
phase [10J. In these magnetic systems it was shown theoretically [11] that in 
the quantum magnetic incommensurably modulated system exist dispersionless 
excitations and that their presence may explain the observed Al-NQR and ther- 
modynamic anomalies observed in CeAl 2 [12] and [15]. and in Ga-NMR study 
of the low-energy excitations in NdGa,2 [H] and [15] . 

In dielectrics the dielectric function behaviour of a composite depends not 
only on properties of infinite and finite clusters, but also on the particles from 
which the composite is formed as well as on their possible domain wall structure 
of the order parameter. The size effect on the ferroelectric phase transition 
in PbTiOs ultrafine particles was experimentally studied in [16], namely the 
size dependence of the transition temperature. The difference in the critical 
temperature for the bulk and for the particle with diameter d is found to be 
T CO o — T c {d) = f,_ d •, , where the constant C is C — 588. 5°C nm and the value 
for the parameter of fitting called critical diameter d c of the particle d c = 12.6 
nm for PbTiOz was found. Note that it is not this critical diameter at which 
the ferroelectric phase in the particle vanishes. Here T coo is the bulk transition 
temperature, T c (d) is the particle transition temperature. Similar effects were 
observed (in earlier works) in KDP [17] and in BaTiOz |18| . 

In pjj] the surface effects on the phase transitions in ferroelectrics are studied 
using a phenomenological theory describing the change of the local spontaneous 
polarisation in the vicinity of the free surface of a ferroelectric thin film. The 
film is kept between metallic electrodes. Depolarising field effects reduce the de- 
viation of this local polarisation from its bulk value as compared to the surface 
effects on phase transitions in other systems. The critical exponents describing 
the behavior of the local polarisation in the vicinity of the Curie- Weiss temper- 
ature T c (d) are the same as the bulk exponents. Only the critical amplitudes 
are changed. 

The Landau theory of phase transitions in thick films is developed in |20J 
for a one-component order parameter. The boundary conditions at the sur- 
faces of a film with thickness L are given by means of an extrapolation length 
(5. The exact expressions are given for the critical temperature and the order 
parameters profile in terms of elliptic functions, and the nature of the phase 
transition is discussed. The dependence of the crystal structure on the particle 
size in BaTiO^ is studied in [21]. The Curie- Weiss temperature of BaTiO^ 
powder has been investigated in the particle size range from 0.1 to 1.0 /jm. The 
transformation from the ferroelectric tetragonal to the paraelectric cubic sym- 
metry occurs at the critical particle size of 0.12 \im at room temperature, and 
the Curie- Weiss temperature drops below the room temperature at the critical 
particle size. In [22] the size driven phase transition in ferroelectric particles is 
studied. The spatial distribution and size dependence of the polarisation is cal- 
culated numerically. Theoretical values for the critical diameter d c in BaTiOj, 
and PbTiOs are compared with experimental values. The surface layer thick- 
ness 6 depends on the particle size. The predicted critical size of the particles 
is smaller than that from the experiment. In [23] the size dependence of the 
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ferroelectric transition of small BaTiO^ particles is studied taking into account 
the effect of depolarisation. The depolarisation energy is reduced in a crystal 
that contains domains of different polarisation. Authors are considering cubic 
particles with alternating domains separated by 180° domain walls. The de- 
polarisation energy and the domain-wall energy contributions are taken into 
account in the Landau-Ginzburg free-energy density. Assuming a hyperbolic 
tangent polarisation profile across the domain wall, the domain- wall energy and 
the domain-wall half thickness can be obtained by minimizing domain free en- 
ergy with respect to this thickness. BaTiO^ is not a perfect insulator, therefore 
a Schottky space charge layer on the surface that shields the interior of the 
crystal from the depolarisation field is considered. The equilibrium polarisation 
P and the domain width D w are found. The results show that the ferroelectric 
transition temperature of small particles can be substantially lower than that of 
the bulk transition temperature as a result of the depolarisation effect. Conse- 
quently, at a temperature below the bulk transition temperature, the dielectric 
constant can peak at a certain cube size L. The results agree with the existing 
experimental observations. According to the authors the theory can also be 
applied to other ferroelectric materials such as KH-iPO^ or PbTiO^. Omit- 
ting the gradient term, |24|, in the free energy expansion, the noncrystalline 
surface layer and depolarisation effect lead to the critical temperature depen- 
dence T c (d) = T coo — j^g where D is connected with the correlation length (: 
D = ( 2 cio | T — T c |, ao is the coefficient from the Landau expansion, d is the 
diameter of the particle. Thus a phenomenological study of the size effect on 
the phase transitions in ferroelectric particles was done. Spontaneous polari- 
sation dependence on the diameter d in spherical particles is studied using the 
Landau phenomenological theory. The spatial distribution of the polarisation 
is obtained numerically. A size-driven phase transition is found. In |25j the 
thickness dependence of the dielectric susceptibility of ferroelectric thin films 
was studied within a phenomenological theory. If the spontaneous polarisation 
is reduced in the surface layer, the mean susceptibility of the film increases with 
the decrease of the film thickness. A size-driven phase transition will take place 
at the critical thickness. If the temperature-driven phase transition of the bulk 
is of the second-order type, the size-driven transition will be accompanied by 
a dielectric divergence. If it is the first-order, a finite dielectric peak will ap- 
pear. If the spontaneous polarisation is enhanced in the surface layer, the mean 
susceptibility of the film decreases with the decrease in the film thickness. No 
size-driven phase transition and hence no dielectric anomaly will occur in this 
case. In [26] the size driven phase transition in nanocrystalline BaTiO^ was 
studied. Nanocrystalline powders of BaTiOs with a narrow size distribution 
were produced. The powder consisted of the crystallite sizes whose dimensions 
ranged from 10 nm to 1 \im. The paraelectric-ferroelectric phase transition was 
found to disappear below a critical crystallite size of 49 nm. In accordance 
with this observation a structural transition towards a cubic symmetry became 
apparent in the x-ray diffractograms. 

The size-induced diffuse phase transition in the nanocrystalline ferroelectric 
PbTiOz has been studied experimentally in [27] ■ The size effects were found 
to become important only below 100 nm (the coherently diffracting x-ray do- 
main size) . The tetragonal distortion of the unit cell related to the spontaneous 
polarisation decreases exponentially with size and vanishes at 7 nm. T c (d) de- 
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creases gradually but the transition becomes increasingly diffuse alongside of a 
diminishing size from 80 to 30 nm. The ferroelectric ordering probably persists 
down to 7 nm. 

The size effects on cells in ferroelectric films are studied in [28] within the 
phenomenological theory. The lateral size dependence of the Curie- Weiss tem- 
perature and the polarisation is obtained and the stability of the ferroelectricity 
of the cells is discussed. The phase-transition behavior of the spontaneous po- 
larisation and susceptibility of the ferroelectric thin films is studied in [29] that 
is based on the Ising model in the transverse field. Modification of the ex- 
change constant and the transverse field in the surface layer might lead to the 
spontaneous polarisation and the Curie- Weiss temperature changing in a differ- 
ent direction as well as in the same direction according [29]. The Curie- Weiss 
temperature of the films can be enhanced, then the temperature dependence 
of the spontaneous polarisation shows a tail-like structure. There exist two 
peaks in the susceptibility-temperature curve in this case. One is located at the 
Curie- Weiss temperature of the films. The intensity of this peak decreases with 
increasing film thickness. The other peak is around the bulk Curie- Weiss point. 
The intensity of this peak increases with increasing film thickness. The theory 
[29J gives a reasonable description of the experimental facts for the triglycine 
sulfate (TGS) thin films. 

Surface effects and size effects for ferroelectrics with the first-order phase 
transition are studied in [30]. The contribution of surface to the free-energy 
expression was studied for the transverse Ising model with four-spin interac- 
tions taken into consideration. It was shown that a P 4 term should be added 
to the surface terms in the free-energy expression. The surface and size effect 
on polarisation and Curie- Weiss temperature were studied using the newly de- 
veloped free-energy expression. Experimental results were discussed using the 
free-energy expression. In [31] the finite size effect in ferroelectrics results in a 
structural instability. This leads to a limitation of physical sizes and dimensions 
of materials in which electric dipoles can be sustained. The size dependence on 
the Curie- Weiss temperature is calculated with consideration of crystallographic 
anisotropy. The mean-field theory gives that the limitation of the transverse 
critical sizes is strongly dependent upon the thickness. The limitation of the 
critical thickness is also closely associated with the transverse critical sizes. 

In [32] a new phenomenological theory of size effects in ultrafine ferroelectric 
particles of PbTiO^ is discussed. This model is taking size effects on the phe- 
nomenological Landau-Ginzburg-Devonshire coefficients into consideration and 
can successfully explain the observed size effects on the Curie- Weiss tempera- 
ture, c/a ratio, and thermal and dielectric properties of the lead-titanate-type 
ferroelectric particles. Theoretical and experimental results for PbTiO^ fine 
particles are compared and discussed. The relationship between this model and 
the model of Zhong et al. [24] is discussed. In [33] the size effects on ferro- 
electricity of ultrafine particles of PbTiOz have been studied by high resolution 
transmission electron microscopy. The diameter of the ultrafine ferroelectric 
lead titanate particles is ranging from 20 to 2000 nm. The crystal structure, 
surface morphology, domain- wall structure and surface reconstruction have been 
studied. All the particles had tetragonal structure. The lattice constants c/a 
ratio and the domain size decrease with decreasing particle size. The particles 
became monodomain when their diameter was 4.2 - 20 nm. A domain wall 
width of 14Awas deduced for the 90° domain walls. 
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Grain-size effects on the ferroelectric behavior of dense nanocrystalline BaTiO^ 
ceramics are studied in [34]. A progressive reduction of the tetragonal distortion, 
the heat of transition, the Curie-Weiss temperature and the relative dielectric 
constant have been studied with grain size decreasing from 1200 to 50 nm. Ac- 
cording to the authors, the correlations between the grain size, the tetragonal 
distortion, and the ferroelectric properties strongly support the existence of an 
intrinsic size effect. From the experimental trends the critical size for disap- 
pearance of ferroelectricity has been evaluated to be from 10 to 30 nm. The 
strong depression of the relative permittivity observed for the nanocrystalline 
ceramics is ascribed by the authors to the combination of the intrinsic size ef- 
fect and to the size-dependent dilution effect of a grain boundary dead layer. 
Size effects on Curie- Weiss temperature of ferroelectric particles were studied in 
|35j within a model without any free adjustable parameters. The model predicts 
that T c decreases with decreasing particle size. The predictions of the model are 
in agreement with experimental results for PbTiO^ and BaTiO^. For similar 
studies see also [36] . 

Influence of the domain wall structure on the dielectric response in KH2PO4 
type crystals was studied in |37J, where one can find also references to domain 
wall dynamics effects in KDP and CSH2ASO4. In this type of the materials the 
domain walls freeze- in below some temperature and a multirelaxation behaviour 
of the dielectric function is observed. The Landau theory of 180° domain walls in 
BaTiOz type ferroelectric particles was studied recently in [38] for the particles 
of the rectangular form. It was assumed that the surface charge is completely 
compensated. The experimentally found temperature dependence on the size 
of the particle is more precisely described by a domain wall if the sixth-order 
term in the Landau free-energy expansion is taken into account [38]. Dielec- 
tric constant and correlation length of ferroelectric particles in the paraelectric 
phase was studied in [39] . The finite size of a particle influences the temperature 
dependence of the dielectric constant and the correlation length. The deviation 
from the Curie- Weiss law in the ferroelectric particles originates from the de- 
crease of the long-range correlation. Glinchuk et al. [ID] studied variationaly 
the depolarisation and the surface tension effects for polarization, susceptibility 
and critical temperature T c {d) diameter dependence for spherical particles. 

In this paper we are studying other sources of the observed complicated mul- 
tirelaxation behaviour of the dynamic susceptibility in ferroelectric-dielectric 
composites by analytical calculations. Phase transition in a spheroidal particle 
from the paraelectric to the ferroelectric phase as well as a dynamic suscepti- 
bility are affected by the spheroidal (more generally by ellipsoidal) shape of the 
ferroelectric particles. Assuming that the surface charge had been compensated, 
the spheroidal shape does not lead to plane polarisation waves due to the sur- 
face effect: the boundary condition in a spheroidal particle with a surface layer 
lead to a multirelaxation response. There is an infinite number of responding 
modes, thus there is also an infinite number of resonance frequencies. It is the 
aim of this paper to discuss namely these effects of the shape of a ferroelectric 
particle on its dielectric response. Firstly we describe a model which enables 
us to calculate a dynamic dielectric function within quasistatic approximation. 
Deriving dynamic equations for the order parameter, the electric polarisation 
vector components, we discuss the paraphase properties of the response. As a 
result we find multirelaxation behaviour of the dielectric susceptibility in the real 
representation. It is then easy to find frequency dependent form of the dielec- 
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trie susceptibility of a spheroidal particle. Its form enables us to find resonance 
frequencies, their number is infinite. Consequently the transition temperature 
dependence on the particle diameter and on other parameters of the model may 
be found easily and compared with the experimental results. 



2 Model 

Let x be a generalized dipole length of the normal polarisation mode with an 
effective charge q. The one-component order parameter P, electric polarisation 
component, has the form P=xq. Let m be an effective mass of the mode. The 
relaxation equation has the form which follows from the free energy F for the 
systems of the fourth order in the order parameter P, see in |22j : 



F = J dv[^P 2 + f P 4 + §(VP) 2 ] + J da[~] - J dv[E.P] (1) 

The first part of the equation |T]) is a volume contribution to the free energy, 
the second one is a surface contribution to the free energy. Last term is due to 
the applied electric field E acting on the electric polarisation P. Here A, B, and 
D are usual parameters of the Landau free energy expansion, E is an electric 
field, 6 is a surface parameter defined by, see in |24j : 

I 5J-4J S 

5 a J { ' 

where J is a ferroelectric interaction constant in the bulk, J s is a ferroelectric 
interaction constant on the surface, a Q is the bulk lattice constant. In (f2|) it is 
assumed for simplicity that the bulk lattice structure is a cubic one. Note that 
the limit 5 = corresponds to strong (anti-) ferroelectric interactions on the 
surface of the particle. On the other hand for strong bulk ferro- (antiferro-) 
interactions J the constant S remains finite and nonzero. For the ferroelectric 
particles with J = J s the constant S is independent on the characteristic size of 
the film. In spherical particles the constant S is dependent on the diameter d of 
the particle: 

where S m is the extrapolation length for d = oo. As noted in [24], to get the 
spatial distribution of the polarisation it is necessary to minimise the free energy 
(HJ) for a given boundary conditions which is difficult to solve analytically. The 
authors of [23] found the spatial distribution of polarisation numerically. We 
have found this distribution analytically for very large values of the ratio -j, 
for which the surface contribution to the free energy (JTJ) is minimised by zero 
surface polarisation. 



3 Equations for the Soft Mode 

The dynamic equation for the generalized coordinate x has the form: 
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m W X + T dl X = ~ q{ Jp -d-rSP ) = qE '" (4) 
where E e ff is an effective electric field acting on the charge q, here P = q.x 
is the polarisation dependent in the position r in the particle. Using the free 
energy F form from ([I]) we find that the equation of motion (J4j) has the form: 

8 2 8 8 2 P 1 BP 

+ T dl X = qE + q( - AP ^ + + r^ )] (5) 

It is convenient to introduce the mass m* related to the charge q and the 
relaxation constant T related to the charge q: 

m 

m = — 
T 

The equation of motion has now the form using ((6]) : 

d 2 d d 2 P 2dP 

ot at or 2 r or 

The form (7]) of the equation of motion is the form which we will use in next 
sections. The boundary condition for a particle with the diameter d has the 
form [2l] : 

dP P 

^ + 7 = ° (8) 

at the surface r = | . It is further assumed that the polarisation P(r) in J7]) 
remains finite at the centre of the particle lim^o P(r) < oo. Note, that the 
equation ((8]) is an algebraic equation relating the P value at the point r = |, 
and the ^ value at the same point. In the limit = 0, the surface boundary 
condition is P(r = |) = 0. This boundary condition will be used in this paper. 
This condition is also compatible with existence of the surface layer with zero 
electric polarization (dead layer) which may exist around the polarized region 
of the particle. 



4 Polarisation Mode in the Paraelectric Phase 

The fourth-order term in the free energy (T]) may be neglected in the mode 
analysis if the temperature T is above the Curie- Weiss temperature T c at which 
the paraelectric phase transforms to the ferroelectric one. In our model above, 
the free energy corresponds to a displacive type of ferroelectrics which undergo 
a second order phase transition described by a one dimensional order parameter 
P. The order-disorder case has to be treated in another way. 

In the free energy F the constant A is of an obvious form A = a\(T — T c ), 
with the positive constant a\ > 0. Let us assume that our particle of the 
spherical shape is under influence of the time dependent spatially homogeneous 
electric field E = E^.e-xjpljuit) oscillating with the frequency to, here j is an 
imaginary unit. Thus we assume that the wavelength of the electric field is 
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much larger than the particle characteristic size. Then electric polarisation 
P = P(r).exp(ju)t) has an amplitude, complex in general, which is dependent 
on the position in the particle, and its time dependence is of the same type as 
of the external electric field. The dynamic equation of motion J7]) has the form 
of a second order differential equation for the polarisation P(r) : 

8 2 P 1 BP 

(-mV + jT*u + A)P(r) =E U + D{— + - — ) (9) 

Substituting P(r) = -^u(r) into the equation ([9]), where u(r) is an unknown 
function, we find the following form of the equation of motion J9]) for u(r) : 



( -m*w a + jT*u> + A)u(z) = EMI + 7^(^ + ~1T< Z ) - ^t<*)) 




z dr z 2 

(10) 

Here r = z| where the variable z, < z < 1, was introduced. The 
corresponding boundary condition has a simple form. At the centre of the 
sphere the polarisation remains finite lim r _,.Q < oo. At the surface we have 
u(z = 1) = 0. To solve the equation of motion (fTO)) let us start from the 
well-known Bessel equation of the form [41J and [42J: 

—^ + -—u{z)-^u{z) = -b 2 u{z) 11 
or- 2 z or z z 

The equation 1|11|1 has as its solution the Bessel functions </±i(A,z), where 
b 2 = A 2 . The root A 2 is defined as that point for which the Bessel function 
vanishes J±i(Ai) = 0. 

The equation of the motion (fTt)|) for the order parameter P(r) in its trans- 
formed form P(r) = -^u(r) has the solution u(r) — Y?i=T a iJ+± (A^z) where 
z is expressed through r as it is given above, which is finite in the centre of 
the sphere r = due to the fact that the relevant Bessel functions J + i(Ajz) 

J + i(A;z) r— 

in the corresponding limit are finite for every i, lim z _ vo — — = y Note 

that one finds lim 2 _»n — h^ 1 — = 00, and thus the Bessel functions J_i (XiZ) do 
not satisfy the condition of the finite solution in the centre of the sphere. It 
is well known, that the roots for the Bessel function J + i (Aj) = are Aj = in. 
Expanding in l|10p also the square root of z into Bessel functions: 

4~z= hJ +h (Kz), (12) 

i—l 

where 6 4 = — -^|£(— 1) J , we find that the expansion coefficients a 4 have the 
form: 



-m*uj 2 + jT*u + A+ ^ 
The space dependence of the electric polarisation P(r) is given by: 
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P(r) = X (r,uj)E u (14) 

concerning the amplitude of the polarisation. Here x( r , w ) is the r-dependent 
susceptibility: 

i=oo Af-i V ;+1 sW^M 

X[ ' J r ^ -m* W 2 + ir % w + A + ^ 

which is of a multirelaxation type. As we see, in <fl5|) there is an infinite 
number of the oscillator-like contributions with the corresponding damping, and 
there exists a phase shift between the mode contributions, which is r-dependent. 
Note, that at the centre of the particle the susceptibility remains finite. 



5 Dynamic and Static Susceptibility 

The dynamic susceptibility x( w ) f° r the particle as a whole is that susceptibility, 
which we obtain from (fT5|) by summing up the contributions over different po- 
sitions within the sphere, relating this quantity to a unit volume of the particle: 

xM = 77 I *( r > ^r 2 dr = £ 2 ^g— (16) 

Here V = is the volume of the particle. We see that the strength of 
every oscillator contribution decreases in the power law, i.e. as i~ 2 , for i (an 
integer) increasing. 

The oscillator resonance frequency w% is given by: 



LO 



2 



1 , , 4ir 2 i 2 D, , x 

m* er 

We see that different resonance frequencies increase as i 2 , and the last term 
in (fl7|) decreases as d 2 with increasing diameter d. The relaxation constant 
remains unchanged. 

Static susceptibility may be easily found from (fl6|) and has the form: 



1 



i—oo 6 



Xstat = j} I X{r, w = 0) Wdr = V S (18) 

The limits of zero frequency and of volume average are interchangeable lim- 
iting processes due to the finite volume of the particle. It can be easily seen 
from lfl8|) that the static susceptibility is no more of the simple Curie- Weiss 
type which would correspond to the first term in the sum. 

The critical transition temperature T c (d) can be found from the lowest fre- 
quency mode, the soft mode, and it is given by: 

T c (d) = T C)00 - ^-R (19) 
aid z 

It increases to its bulk value T Ci00 quadratically with the increase of the parti- 
cle diameter d. There is a size driven phase transition with a decreasing diameter 
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d. For the diameter lower than the critical diameter d < d c = 2ir</ — ^ — the 
ferroelectric phase vanishes. 

6 Comparison with Experiments 

We have compared the theoretical dependence lfT9|) found in this paper for the 
critical temperature dependence on the diameter d of the particle with exper- 
imental dependence of the critical temperature on the particle diameter for 
PbTiOz and BaTiOs particles. Experimental results for PbTiO^ are found 
from the measurements in [16]. The authors of [16] used a theoretical formula 
T coo — T c (d) — ^ d 2 d j , where they have found the constant C, C — 588. 5°C nm, 
and the value for the critical diameter d of the particle d c = 12.6 nm to fit the 
experimental points. Thus they have found that there exists a critical diame- 
ter below which there is no ferroelectricity. Experimental results for BaTiOz 
are based on the measurements in [18] and used in (24]. They were described 
by the theoretical formula T c (d) = T coo — |^ where D' is connected with the 
correlation length £: Z> = £ 2 ai | T — T c \, ai is the coefficient from the Landau 
expansion, d is the diameter of particle. The critical diameter was found to be 
44 nm theoretically and 115 nm experimentally for BaTiOz- The same authors 
have found that for PbTiOz the value is 4.2 nm (theoretical value) and 13.8 nm 
(experimental value from [16]) and 9.1 nm (experimental value from [24]). It is 
interesting that Glinchuk et al. |l0] obtained a formula for the critical temper- 
ature T c (d) dependence on the diameter of a particle by a variational method 
which has the form T c (d) = T CjOC (l — % — jz^t) where dj, is the correlation 
diameter renormalized by the surface tension, d c is the correlation diameter 
renormalized by the surface tension and by the depolarisation field and d\ is a 
parameter. As we can see, the surface tensions and depolarisation fields lead to 
the dependences of the critical temperature used in the papers above as theoret- 
ical formulas for fitting the experimental points. Our calculations in this paper 
of the critical temperature dependence on the particle diameter are exact. It 
is known that variational calculations always lead to a higher free energy than 
exact calculations, thus our formula (fl9|) (depolarisation field is assumed to be 
zero) which leads to dependence instead the Glinchuk et al ^ dependence, is 
expected to be better for neglected depolarisation effects. If the depolarisation 
effects are taken into account we may expect that a constant d m in our formula 
(jX9j) will be introduced in ^ d _ l d p instead of as a similar dependence is 
found in the Glinchuk et al. formula. This leads to the following fitting formula 
for the critical transition temperature T c (d): 

r c (d) = r C|0O (i- d \ ) (20) 
(rf- d m y 

This formula is expected to correspond to the Glinchuk et al. formula [40J 
in some sense. In this formula, the depolarisation and surface tension effects are 
taken into account for diameters near the corresponding correlation diameter. 
Our above-given results were obtained using the exact calculation of the polar- 
isation profile induced by an external electric field in the paraelectric phase for 
the large -j limit. The dependences from l|20p with ^ d _ d p and from lfT9|) jp 
are fitted to the experimental data for PbTiO^ , see Figure ([1]), and for BaTiOz 
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, see Figure |(2]). Far the best fit - the curve (1), is such in which we are tak- 
ing the critical temperature T c oo , the parameter d and the parameter d m as 
fitting parameters according to ([20| . We have fitted this dependence (f20|) also 
using the critical temperature as a fixed parameter taking the values from the 
experiment and fitting the parameter d and the parameter d m . We have found 
that such a fit is still better than the fit using the formula (fl9|) . see the curves 
(2) and (3) respectively in Figure ([1]) and Figure |(2|). 

Critical temperature T c in degrees Celsius (we are using degrees Celsius ac- 
cording to available experimental points) dependence on diameter d of particles 
for PbTiOy, is in Figure (1). Theoretical curves are: curve (1) for T C!OQ , d D 
and d m as fitting parameters (we have found T Ci00 = 511. 5°C, d a = 4.8 nm 
and d m — 7.3 nm with quality of fitting procedure given by R 2 = 0.811 and 
d c = 12.1 nm), curve (2) for T c oo and d as fitting parameters with d rn = (we 
have found T Ci00 = 519.0° C and d a = 7.5 nm with R 2 = 0.8024 and d c = 7.5 
nm), curve (3) for d as a fitting parameter and d m = (we have found d = 5.9 
nm with R 2 = 0.6425 and d c = d Q ) with T Ci00 fixed to 500. 0°C. In [16] regard- 
ing PbTiOz ceramics the diameter corresponding to our d m diameter is found to 
be 12.6 nm, compared to our dead layer width 7.3 nm in Figure (jTJ) . The critical 
diameter corresponding to our d c was not evaluated in [16J. Note that the dead 
layer diameter (=2 x dead layer radius) is in JT6] interpreted as a critical size 
below which ferroelectricity becomes unstable. 

Critical temperature T c in Kelvins (according to available experimental points 
we are using Kelvins) dependence on the diameter of particles - BaTiO^ . Theo- 
retical curves are: curve (1) for T c oo , d Q and d m as fitting parameters (we have 
found T Cj00 = 390.1 K, d = 12.1 nm and d m = 92.8 nm with a quality of the 
fitting procedure given by R 2 — 0.9938 and with d c — 104.9 nm), curve (2) for 
T Cj00 and d Q as fitting parameters with d m — (we have found T Cj00 = 410.2 
K and d a = 56.5 nm with R 2 = 0.7768 and d c — d a ), curve (3) for d Q as a 
fitting parameter with fixed T c>00 = 391.0 K from [24] and d m = nm (we have 
found d = 49.4 nm with R 2 = 0.6940 and d c = d Q ). In [H] the theoretically 
predicted value of the critical diameter is 44 nm and experimentally found 115 
nm, they should be compared to our 104.9 nm for the curve (1), 56.5 nm for 
the curve (2) and 49.4 nm for the curve (3). We see that our value of the 
critical diameter from the curve (1) is near the experimental value. Our values 
for curves (2) and (3) are near the theoretical value (for which there no dead 
layer was assumed in [21] )■ In [HI] on dense BaTiOz ceramics with grain size 
decreasing from 1200 to 50 nm, the authors have found strong support for the 
existence of an intrinsic size effect. From the experimental trends the critical 
size for disappearance of ferroelectricity has been evaluated to be from 10 nm 
to 30 nm. The strong depression of the relative permittivity observed for the 
nanocrystalline ceramics can be ascribed to the combination of the intrinsic size 
effect and of the size-dependent dilution effect of a grain boundary dead layer. 
We have found that the quantity d m has the value 92.8 nm for the curve (1). 
This quantity, which contributes to the total diameter of the particle, leads to 
the part of the diameter which corresponds to polar region with d=11.9 nm for 
the critical diameter d c — 104.9 nm which we have found, see Figure ([2]). Thus, 
if we identify d m with the dead layer diameter, then we see that the critical 
diameter found by our fitting procedure is 3 to 10 times larger than that found 
experimentally. 
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7 Conclusions 



A phase transition in a spheroidal particle from paraelectric to ferroelectric 
phase as well as dynamic susceptibility are affected by the spheroidal shape 
of ferroelectric particles. This shape leads to many resonance frequences in the 
dynamic susceptibility. Sources of the observed complicated multirelaxation be- 
haviour of the dynamic susceptibility in ferroelectric-dielectric composites are 
of a different type, the shape of particles in microcomposites is one of them. In 
our model in this paper we assume that the surface charges are compensated. In 
general, there is a different polarisation at the surface of the particle from than 
in the bulk. The spheroidal shape does not lead to plane polarisation waves 
due to the surface effect: the boundary conditions in a spheroidal particle with 
a surface layer lead to a multirelaxation response. We have found that for the 
surface boundary condition P = on the surface there is an infinite number 
of the responding modes as well as infinite number of resonance frequencies. 
We have found formulas which enable to discuss the effects of the ferroelectric 
particle shape (spheroid) on its dielectric response. First we described a model 
from which we have found the dynamic dielectric function within the quasistatic 
approximation. Dynamic equations for the order parameter, the electric polar- 
isation vector component, were found. We discussed the paraphase properties 
of the dielectric response of the particle. The Cole-Cole diagram is not of the 
semicircle form. The static susceptibility within the paraphase is discussed. It 
does not have the Curie- Weiss like behaviour as in the bulk. 

The transition temperature decreases with the decreasing diameter d of the 
particle quadratically. We discussed our theoretical results comparing them 
with experimental results for the critical temperature diameter dependence in 
PbTi0 3 and in BaTi0 3 . We have found that the formula T c (d) = T Ci00 (l - 

d 2 

(d-d ) 2 ) wnere ^c,oo is a fitting parameter in general together with d a and d m , 
the last parameter is a parameter describing the dead layer thickness. 
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Figure 1: Critical temperature T c in degrees Celsius (according to experimental 
points - full triangles) dependence on diameter of particles - PbTiO^. Theoret- 
ical curves are: (1) for T c >00 , d and d m as fitting parameters, (2) for T c oo and 
d Q as fitting parameters with d m = nm, (3) for d as a fitting parameter with 
Tcoo = 500°C and d m = nm. 
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Figure 2: Critical temperature T c in Kelvins (according to experimental points - 
full traingles) dependence on diameter of particles - BaTiO^. Theoretical curves 
are: (1) for T Cj00 , d and d m as fitting parameters, (2) for T c ,oo and d as fitting 
parameters with d m = nm, (3) for d as a fitting parameter with T C]00 = 391 
K and d m = nm. 
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